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Abstract 

The most fundamental characteristics of a physical system can often be deduced from 
its behaviour under discrete symmetry transformations such as time reversal, parity 
and chirality. Here we review basic symmetry properties of the relativistic quantum 
theories for free electrons in (2 + 1)- and (l + l)-dimensional spacetime. Additional 
flavour degrees of freedom are necessary to properly define symmetry operations in 
(2+1) dimensions and are generally present in physical realisations of such systems, 
e.g., in single sheets of graphite. We find that there exist two possibilities for defining 
any flavour-coupling discrete symmetry operation of the two-flavour (2+ l)-dimensional 
Dirac theory. Physical implications of this duplicity are discussed. 
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1. Introduction 

The realisation of a two-dimensional (2D) quasi-relativistic electron system in single- 
layer graphene [ 1 ] and at the surface of topological insulators [2] has renewed interest 
in 2D versions of the Dirac equation [3, 4] for a free particle, 

H ifr = id,ifr, (la) 

with Hamiltonian 

# D (p) = v-p+/?m. (lb) 

We summarise here some features arising due to the reduced dimensionality, 
focusing especially on symmetry properties. Comparison is made between previously 
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considered (2+l)-dimensional quantum electrodynamics (QED3) [5-7] and quasi- 
relativistic condensed-matter systems. [8-13] 

1.1. Dirac Equation: Basics We employ the general expression for the Dirac 
Hamiltonian given in Eq. (lb) in terms of the velocity operator v and the matrix valid 
in any spatial dimensions d > 1 . Our notation emphasizes the fact that the matrices a 
(= v) commonly used for the Dirac equation in (3 + l)D have the physical interpretation 
of being the operator for the velocity of a free particle. While the explicit form of these 
matrix operators depends on the dimension, the structure of the Hamiltonian Hd does 
not. As the Dirac Hamiltonian for a free particle commutes with momentum p, we 
can choose plane waves as basis in real space and represent p as an ordinary vector in 
R d . In the following, we adhere to this convention of condensed-matter physics, which 
implies that only the constant matrices v and are quantum-mechanical operators. We 
use units where the Planck constant fi = 1 and also |v| 2 = 1. 

The spectrum of the Dirac Hamiltonian (lb) has two branches, one at positive and 
one at negative energies. These are given by 

£ ± (P) = ±VlPl 2 + m2 (2) 

independent of the number d of spatial dimensions. The eigenstates in position 
representation have the general form 

^(r) = e'P r 4 . (3) 

In (3+l)D, the x% are four-component spinors that encode the energy -band index 
and the Dirac particle's intrinsic angular-momentum (spin) vector, where the 
corresponding spin operator is defined as 

S--vxv. (4) 
Ai 

Here the symbol 'x' denotes the cross product, defined formally as for vectors in R 3 . 

The relativistically covariant form of the Dirac equation is obtained by multiplying 
Eq. (la) by f5 and re-arranging to exhibit derivatives w.r.t. time and spatial coordinates 
in a unified fashion. Using the definitions y° = /?, y ; = /? vj with j = 1, . . . , d, as well 
as the components = id^ of the (d + 1)- vector of momentum [14], one finds 

(y Vn -mjifr = , (5) 
where we have used the convention that repeated indeces are summed over. 

1.2. Symmetries The symmetries of the Dirac equation include operations that 
reflect the properties of Minkowski space time. (Point group operations give rise to 
the Lorentz group, whereas the Poincare group also includes translations in space and 
time.) These symmetries are also relevant for a classical (non-quantum) description 
of relativistic systems. On the other hand, in general we also have symmetries that 
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are present only in a quantum mechanical description of relativistic systems. Each of 
these categories include discrete and continuous symmetries. 

In the following, we will focus mostly on the discrete symmetries that are specific 
for the quantum-mechanical Dirac equation. Among the spatial symmetries, we will 
concentrate on the nontrivial operations of parity (space inversion) and rotations by 
In. Moreover, we will consider the operation of time reversal (also called reversal of 
the motion). [15] The Hilbert space of the (single -particle) Dirac equation includes 
solutions with positive and negative energies corresponding to particles and holes (i.e., 
anti-particles). These are related by particle-hole conjugation and energy-reflection 
symmetry. [9, 14, 16-18] Lorentz invariance will not feature at all. 

We continue this Section with a brief introduction to particular discrete symmetries 
and review their representation for the (3 + l)D Dirac case. 

1.2. 1. Spatial Symmetries — in particular, Parity A spatial symmetry transformation 
g is represented by a unitary operator P(g) acting in the Hilbert space of the 
Hamiltonian H, where the set {g} of symmetry transformations forms a group G. To 
describe the symmetry of an observable O under the symmetry transformation g we use 
the short-hand notation gO = P{g) O P~ l (g). When P(g) acts on the basis functions 
{| v)} of H, we can expand P{g) \ v) in terms of the basis functions {| v)}, 

Pig) |v> = 2d(sVI,u> - < 6 > 

so that the set of matrices {T)(g)) forms a representation for G. Then the invariance of 
the Hamiltonian H under the symmetry transformation g implies [19] 

D(g)H(p)D-\g) = H(gi>). (7) 

This equation expresses the fact that we can view the effect of g from two equivalent 
perspectives: we may regard it as a unitary transformation of the basis functions {| v)} 
of H ("active" view). This is equivalent to applying g to momentum p (and position 
r, "passive" view). If g is a symmetry of H, both views must yield the same result. 
The invariance condition (7) holds generally for continuous symmetry transformations 
(such as rotations) as well as discrete operations (e.g., parity). 

Equation (7) includes as a special case parity P (i.e., space inversion). In d = 3 
spatial dimensions, it is characterized via the relations 

9 r = -r, V p = -p (8a) 

TKP) H(p) £T 1 (p) = H(P p) = //(-p) (8b) 

For the Dirac Hamiltonian in (3+l)D, the representation matrix for the parity operation 
is 

DCP)=(3. (9) 
As to be expected, we have T){P 2 ) = +1. 
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7.2.2. Rotations K In (3+l)D, the Hamiltonian H D commutes with the operator of 
total angular momentum J = L + S, where L = r x p is the operator of orbital angular 
momentum, and the spin operator S is given by Eq. (4). Total angular momentum J is 
the generator for finite rotations. Rotations K„ by 2n about an arbitrary axis n maps 
the system onto itself with a minus sign 

£>C%(2tt)) = exp(27n n • J) = -1 . (10) 

In that sense, the symmetry group becomes a double group which is common for half- 
integer spin systems [20]. 

1.2.3. Time Reversal Time reversal (TR) is an anti-unitary transformation and can 
be represented by the operator 

e = Tn< (11) 

where K denotes complex conjugation, and the operator T is the unitary 
transformation that relates the time -reversed (i.e., complex-conjugated) basis spinors 
to the original basis [19, 21]. TR invariance requires that the Dirac equation has the 
same solution as the equation obtained by applying 6 to both the basis functions and 
all tensors that represent physical quantities (such as momentum p), giving 

D(T) H*(p) D(T)- 1 = H{J9 p) = ff(-p). (12) 

Here, * denotes complex conjugation. The minus sign reflects the fact that momentum 
p is odd under reversal of motion. This implies that 6 reverts the velocity v, 
D{T) v* T){T)~ l = -v, but does not change the mass term, £>(T) p* D(T)~ l = /?. For 
the Dirac Hamiltonian in (3+l)D in the standard representation [14] where only v 2 is 
imaginary, we have 

TKT) = -i Vl v 3 ^yV , (13) 

which implies D(6 2 ) = -t. The remarkable property that the representation matrix 
for the time reversal operation squares to the negative of the unit matrix for a (3 + l)D 
Dirac particle arises as a consequence of its intrinsic half-integer spin [22]. 

1.2.4. Particle-Hole Conjugation Particle-hole conjugation is also an anti-unitary 
operation [14, 16-18]. We can express it using the operator 

k = CK , (14) 

with a unitary operator C. Particle-hole symmetry of a Hamiltonian holds if the relation 

£>(C) ff*(~P) 0{CT x = -H(P) (15) 

is satisfied. The operation of particle-hole conjugation reverses momentum and 
spin but leaves position invariant [14]. For the Dirac Hamiltonian (lb) in p 
representation, the relation (15) is equivalent to the conditions D(C) v* T){C)~ X = v 
and D(QP* D(C)" 1 = -yS. For the Dirac Hamiltonian in (3 + l)D, we have [14] 

D(C) = i/3v 2 = iy 2 , (16) 

which implies D(k 2 ) = +1. 
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1.2.5. Chirality and Energy Reflection Symmetry In (3+l)D, an additional matrix 
y 5 = j'y°y 1 y 2 y 3 can be defined that anticommutes with all other y. In the massless 
limit m = 0, y 5 commutes with the Dirac Hamiltonian 

[// D (p,m = 0),y 5 ]=0, (17) 

so that it becomes the generator ^ for a continuous chiral symmetry [14] with T)(x) — 
J 5 - 

In the context of random-matrix theory (RMT), the term chiral symmetry has 
departed from its original meaning in relativistic quantum mechanics [17, 18]. The 
chirality matrix y 5 anticommutes also with the Dirac operator D = y M - m, see Eq. 
(5). Except for zero modes, the eigenfunctions of D occur in pairs u„, y 5 u n with 
opposite eigenvalues (A„, -A n ). This notion [23] has been transfered to other systems. 
In RMT, a chiral symmetry [17, 18] of a Hamiltonian H is associated with a unitary 
operator M that satisfies 

D(M)H(p)!D- l (M) = -H(p) . (18) 

This operation reverses the energy. In recent studies [9], the operator M has therefore 
been called energy-reflection (ER) symmetry. In order to distinguish M from the 
chiral symmetry (17), we will follow the latter naming convention. The operator M 
reverses the velocity, while leaving momentum and position invariant: M p = p and 
Mr = r. For the Dirac Hamiltonian, the corresponding representation matrix D(M) 
must therefore anticommute with the matrices v and (5. 

For the (unperturbed) Dirac equation (1), time reversal symmetry (12), particle- 
hole conjugation (15) and the energy reflection symmetry (18) are not independent 
symmetries. [18] We have 

M = k6 (19) 

giving 

£>(M) = £>(C) [D(T)T , (20) 

and explicitly in (3+l)D 

D(M) = -iy^^-ipDfx) , (21) 
which implies D(M 2 ) = +1. 

2. (2+l)-dimensional Dirac theory 

In this Section, issues relating to the discrete symmetries for the (2+l)D Dirac 
theory are discussed. In (2+l)D spacetime, the Dirac Hamiltonian becomes a 2 x 2 
matrix. In this case, we have only four linearly independent basis matrices that we may 
choose as Pauli matrices [22] o~ x , cr y , and cr-, as well as ctq = 12x2- Yet we still need 
to represent intrinsic angular momentum (spin) and two energy bands. The reduced 
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number of basis matrices thus gives rise to pathologies that prevent the consistent 
definition of many symmetries, except in the mass-less limit (m = 0). 
Explicitly, we can choose a representation where 

v = (cr r , cr y ) and p = cr- . (22) 

Our representation is equivalent, but not identical, to previously discussed ones [5, 6, 
8], as a look at the y matrices shows: 

y° = cr- , y 1 = icr y , y 2 = -icr x . (23) 

2.1. Angular Momentum and Rotations P. Straightforward calculation establishes 
that the (2+l)D version of commutes with the operator of total angular momentum 
J- = L- + 5 - , where L z = x p y - y p x is the z component of orbital angular momentum, 
and 

S, = I (vxv)-z = i(r : . (24) 

Thus the intrinsic spin degree of freedom formally emerges in complete analogy to the 
(3+l)D case; see Eq. (4) and Ref. [4], but only its z component is a relevant quantity 
in (2+l)D. 

Rotations % z by 2n about the z axis map the system onto itself with a minus sign 

D(R z (2n)) = exp(2m7,) = -1 , (25) 

similar to Eq. (10), so that once again the symmetry group becomes a double 
group [20]. 

2.2. Parity P In three spatial dimensions, parity V is defined as space inversion 
r — > -r. Its interesting aspects arise from the fact that the determinant of the associated 
representation matrix is det(£)(!P)) = — 1, so that parity cannot be expressed in terms 
of infinitesimal unitary transformations. In two spatial dimensions, the transformation 
r = (x, y) — > — r = (— x, -y) is not interesting because it is equivalent to a rotation by n 
and the determinant of the representation matrix is +1. This is why parity in 2D is 
usually defined [5-8] as either 

r = (x, y) -> P x r = (je, -y) (26a) 

or 

r = (x,y)^P y r = (-x,y) . (26b) 

The same transformation laws apply for linear momentum p and velocity v. Clearly, 
parity invariance as a property of a physical operator should be established with respect 
to both operations (26a) and (26b) to count as a proper symmetry. The above 2D 
versions of parity imply that L-, is odd under parity, 



P x L z — Py L z — L- 



(27) 
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If we require that v ■ p is even under parity, we find 

D(P v ) = o- v with v = x,y . (28) 

Furthermore, this implies [5] that S z and J. are likewise odd under parity. In contrast, 
yS must be even under parity. As S z = ft/2 for the (2+l)D Dirac theory, it is impossible 
to satisfy both conditions, i.e., the mass term violates parity symmetry in (2+l)D. Only 
in the mass-less limit m = 0, is invariant under parity. 

We can study parity also for the (/'p(r) given in Eq. (3), which are eigenfunctions 
for both the Dirac Hamiltonian and linear momentum p = (p x , p y ). (See also Problem 
3.4 in Ref. [24].) In the m = limit and for E 2 = p 2 + pi + 0, we have (apart from a 
normalization constant) 

^(r)oce*- r ( £ * ) (29) 



pv \Px + ip 

If parity is defined as above it follows that 

D{P V ) (A*(r) = i% tfr^CP; 1 r) . (30) 

with a phase factor r/p. (Of course, we have P~ x = r Py) Similar to Eq. (7), this 
equation illustrates the fact that we can view a symmetry transformation such as parity 
from two equivalent perspectives: we may regard it as a unitary transformation of the 
expansion coefficients (i.e., the spinors). This is equivalent to the inverse coordinate 
transformation, where position r and momentum p are mapped on the inversely 
transformed quantities. As expected, parity flips position vectors and momenta, while 
the energy is preserved. Also, the expectation value of S z is reversed under parity. 

2.3. Time reversal The basic requirements for TR invariance imply that (3 should 
be even under TR while 5 - should be odd. As in the case of parity, the fact that both 
quantities are represented by o% in the (2+l)D situation makes it impossible to find 
such a TR transformation, i.e., the mass term violates both parity and TR symmetry. 
However, TR can be defined consistently in the massless limit: the usual [22] 

Q = T<K with F>(T) = io- y , (31) 

where o~ v is the imaginary Pauli matrix, implies 

for arbitrary spinors I Gl ] and O{0 2 ) = — 1. More explicitly, we obtain for the wave 
\a 2 ) 

function (29) 

e^=D(T)^ = ^^_ p (33) 

with a phase factor i.e., momentum is reversed, while the energy is preserved. 
Also, the expectation value of S z is reversed under TR. 
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2.4. Particle-hole conjugation C As v x and /? are real matrices in our adopted 
representation of the (2+l)D Dirac theory while v y is imaginary, the matrix C 
entering the definition (14) of particle-hole conjugation needs to commute with v x 
and anticommute with both v y and f3. These conditions are satisfied by the matrix cr x . 
Hence, the particle-hole-conjugation operation k = CK can be consistently defined for 
the (2+l)D case via the representation matrix 

£>(C) = a- x (34) 

so that D(k 2 ) = +1. 

2.5. Energy Reflection Symmetry Ai ER symmetry requires D(M) to anticom- 
mute with all Dirac matrices (i.e., the components of v and fi). No such matrix exists in 
(2+l)D and, hence, ER symmetry cannot be established for the Dirac model. However, 
D(M) = <r z satisfies the ER condition (18) in the mass-less limit with D(M 2 ) = +1. 

2.6. Chiral symmetry x The fact that no matrix anticommutes with all three Pauli 
matrices, i.e., the y^ given in Eq. (23), also implies that no equivalent of y 5 exists 
in (2+l)D [5, 6]. As a result, no chiral symmetry x can he established, even in the 
mass-less limit. 

2.7. Discussion It appears that particle-hole conjugation and spatial rotations are the 
only symmetries that can reasonably be considered within the (2+l)D Dirac theory. 
The fact that consistent parity, TR and ER operations exist only in the zero-mass limit, 
and that chiral symmetry is altogether impossible to define, makes this theory highly 
pathological from a physical point of view. As it turns out, this unsatisfactory situation 
can be remedied by introducing an additional flavour degree of freedom for the Dirac 
fermions [5]. This generalization of (2+l)D Dirac theory will be explored in the next 
Section. 



3. Two-flavour (2+l)D Dirac theory 

It is possible to construct a theory of Dirac fermions in (2 + l)D that exhibits the 
familiar symmetries of parity and TR, even for a finite mass [5]. This theory involves 
four-spinor wave functions describing two flavours (2Fs) of such fermions. The Dirac 
operator (lb) is block-diagonal, and a possible choice for the matrix operators v and /? 
is 

cr x \ / cr v \| , _ /(TV 



"io „,Ho' 'A and -*.)■ <35> 

Such a flavour-symmetric representation was adopted, e.g., in Refs. [25-28] to 
describe the electronic degrees of freedom near the two valleys in graphene. Other 
formulations found in the literature [5-9, 29] are unitarily equivalent to the one used 
here. According to Eq. (24), the spin operator S : emerges as 



1(CT Z 

2 <t z 



Sz = -\ n " ■ (36) 



[9] 
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Thus, in contrast to the single-flavour case, (3 and S z are not proportional to each other 
in the 2F (2+l)D Dirac theory. For completeness, we also give the y matrices: 



KTy 






-io~ v 



-KT X 






iCTy, 



(37) 



3.1. Symmetries Rotations are again generated by J z = L. + S -, which is a 
conserved quantity. Hence, as in the single-flavour case, a 2n rotation is equivalent 
to -1, which is a signature of half-integer spin. 

We define parity, TR and ER as symmetry transformations that connect the flavour 
subspaces. For each of the parity transformations P x ,y, two inequivalent representation 
matrices D ± CP V ) can be realised, 





+cr v 







(38) 



which are distinguished by the property !D ± (P 2 ,) = ±1. Similarly, the anti-unitary 
operator 6 = C T ( K with representation matrices 



£>±<T) 



o 

±iO~y 



-KTy 




(39) 



keeps the 2F (2+l)D Dirac Hamiltonian invariant while reversing velocity, momentum 
and spin. Action on a general state yields 







(-<} 












±a* 2 









(40) 



which implies D±(9 2 ) = ±1. Thus it is possible [30] that Dirac particles in 2F (2+l)D 
behave under TR like spinless particles [having D^O 2 ) = +1] as an alternative to the 
expected spin- 1/2 -behaviour [where = -1]. 

The conditions associated with invariance under particle-hole conjugation k = CK 
are satisfied by 

n ± ( C) ^ 1}) . (4i) 

Unlike the TR and chirality operations, particle -hole conjugation preserves the flavour 
degree of freedom so that we get D±(/e 2 ) = +1. 

Finally, there exist two inequivalent representation matrices for ER symmetry 



D ± (M) = 





+0% 



0~ z 





(42) 



They both anticommute with all Dirac matrices in the 2F (2+l)D model, satisfy the 
relation D E (M) = £>±(C) £> ±E (T)* [s = ±1; cf. also Eq. (20)] and, similarly to the case 
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Table 1 . Summary of 2F (2+l)D operators given in compact notation where Pauli matrices cr, and r, act 
in Dirac and flavour space, respectively, and <xo = To = 12x2- 



V x 


Vy 




7 1 


7 2 


25, 


cr x <8> t 


CT y ® T 




/<x v (gi r. 


-icr x ® t z 


cr_ ® r 


£> + CPv) 








£> + (C) 


»-(C) 


<Ty <8> T x 


iCT v ® T v 




—l CTy <8> 


cr x ® t z 


0-^®To 


£> + (M) 


£>_(M) 










cr z ® t x 


icr- ® r Y 




/cr (gl T y 







of the other flavour-coupling discrete symmetry transformations, are distinguishable 
by the sign of their squares: D±(M 2 ) = ±1. 

A continuous chiral symmetry also exists in 2F (2+l)D models [6, 9]. In our 
representation, the matrices 

«W = U 0°) <«> 

anticommute with ft and commute with v XJ „ thus satisfying the condition for the 
generator of a chiral symmetry. Again, the two matrices are distinguishable by the 
sign of their squares: D±(x 2 ) = ±1. It is possible to express the T)±(x) in terms of 
other Dirac operators and the ER transformations, 

D ± (x) = -iv x v y D ± (M)=l3D*(M) , (44) 

which mirrors Eq. (21). 

3.2. Discussion For reference and to facilitate easier comparison with the literature, 
we provide expressions of relevant operators for the 2F (2 + l)D Dirac model using a 
compact notation (with Pauli matrices crj and tj acting in Dirac and flavour space, 
respectively) in Table 1 . 

The existence of two possible realisations for each of the symmetries S = 
"P v , T, C, M, and x seems puzzling, especially for many of those transformations 
where the squares of the two realisations have opposite signs. However, for 
any specific physical realisation of a 2F (2+l)D Dirac system, the transformation 
properties of the basis functions will uniquely determine for each S which of the two 
operators T)+(S) corresponds to the actual symmetry in this system. For fermionic 
(spin-1/2) particles, TR will be represented by [22] D-(T~). In contrast, D+(l~) applies 
to electrons in graphene [27] whose low-energy band structure (neglecting the real 
spin) realises a 2F (2+l)D Dirac model, with the valley index being the flavour degree 
of freedom [1]. 

In previous discussions [6] of chiral symmetry in 2F (2+l)D Dirac systems, the 
existence of two possible representations T)±(x) was seen to imply the existence of 
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a U(2) symmetry with generators {1, D-(x), D+(x)D-(x)}- Considering also 

the other pairs of representation matrices !D+(S), we find more generally 

D + (S)D-(S) = [~q ° t J . (45) 

for each symmetry S. Hence the product of the two possible representation matrices 
for each symmetry acts differently in the subspaces associated with the two flavours. 
Therefore, in situations where the two flavours are physically indistinguishable (which 
is the case, e.g., in an ideal sheet of graphene), the two representation matrices £)±(S) 
cannot be interpreted as distinct discrete symmetries of this system. 

4. (l+l)-dimensional Dirac theory 

For comparison, we briefly discuss the (1+1)D Dirac model where p = p in H D 
given by Eq. (lb). For definiteness we choose v = cr x and = cr.. Parity requires 

£KP) = cr. (46) 

with DCP 2 ) = +1. For TR, we find 9 = T K with 

£>(T) = cr, . (47) 

Thus we have D(8 2 ) = +1, consistent with the fact that a Dirac particle in (1 + 1)D is a 
spinless object [4]. Particle-hole conjugation is realised by k = CK where 

D(C) = cr x (48) 

with D(k 2 ) = +1. The ER transformation with 

D(M) = -i o-y (49) 

and D(M 2 ) = +1 satisfies the condition (18), and D(jc) = cr x defines a chiral 
symmetry. 

5. Conclusions & Outlook 

We investigated basic properties of discrete symmetry operations for low- 
dimensional Dirac models. An additional flavour degree of freedom is needed to 
introduce the usual discrete symmetries known from Dirac particles in three spatial 
dimensions. For the two-flavour (2+l)-dimensional Dirac model two inequivalent 
representation matrices exist for each symmetry operation. With the exception of 
particle-hole conjugation, the two representation matrices for each symmetry are 
distinguishable by the sign of their squares. 

Symmetry analyses such as that presented in this work enable the classification 
of perturbations to the (2+l)-dimensional Dirac model [18], making it possible to 
draw conclusions about how the perturbations affect the spectrum and, ultimately, the 
physical properties of condensed-matter realisations such as graphene sheets. The 
origin of such perturbations could be disorder [25, 31, 32] or p-dependent corrections 
to the Hamiltonian reflecting the lower crystal symmetry of, e.g., the graphene 
lattice [21]. 
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